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Abstract. We have modeled the photochemical evolu- 
tion of the sulfur-containing species that were observed 
in Jupiter’s stratosphere after the SL9 impacts. We find 
that most of the sulfur is converted to Sg in the first few 
days. Other important sulfur reservoirs are CS, whose 
abundance increases markedly with time, and possibly 
H 2 CS, HNCS, and NS, whose abundances depend on 
kinetic reaction rates that are unknown at the present . 
We discuss the temporal variation of the major sulfur 
compounds, make abundance and compositional pre- 
dictions useful for comparisons with observations, and 
discuss the possible condensation of sulfur-containing 
species. 

Introduction 

Observational reports indicate that many new molec- 
ular species were deposited in Jupiter’s stratosphere 
following the Comet P/Shoemaker-Levy 9 (SL9) im- 
pacts. These new species derive from highly processed 
cometary material and tropospheric Jovian air that 
were modified during the impact explosion, fireball, and 
re-entry shock (‘'splashdown”) of the plume back into 
the upper atmosphere. To examine the photochemi- 
cal evolution of the species that were observed or pre- 
dicted to have formed following the plume splashdown 
(~ 1 0- 1 5 minutes after the actual impact), we use a one- 
dimensional (1-D) chemical kinetics and diffusion model 
similar to the one described in Yung et al. [1984]. Our 
results regarding the sulfur-containing species are pre- 
sented at this time — results regarding nitrogen and 
oxygen species are presented in a companion paper 
[Moses et a/., 1995]. 

The Hubble Space Telescope (HST) Faint Object 
Spectrograph (FOS) obtained ultraviolet spectra within 
the central dark core of the G (and L) site over a period 
of several days and weeks following the impacts [Noll 
et a/., 1995; Ye lie and McGrath , 1995]. These high- 
quality spectra indicate the presence of several sulfur 
compounds at the site; S 2 , CS 2 , CS, and H 2 S are all 
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identified; notably absent are signatures for SO 2 and 
SO. Although SO 2 has been detected in inid-infrared 
spectra [A. Sprague , personal communication, 1994], its 
inferred abundance is several orders of magnitude be- 
low that of the other observed sulfur species. These 
observations suggest that oxygen compounds may not 
have played a large role in the impact chemistry [cf. 
Zahnle et al. 1995]. On the other hand, several Earth- 
based observations imply that oxygen species such as 
H 2 O and CO were abundant at the same G impact site 
[Lellouch et al., 1995; Bjoraker et al., 1994; Meadows 
et al., 1994]. Because the reason for the discrepancy 
between the HST and other observations is not known, 
we concentrate on the dark core regions observed by 
the FOS and assume that the oxygen compounds re- 
main within the HST upper limits. 

The complex nature of the impacts makes it impossi- 
ble for 1-D photochemical models to reproduce all the 
phenomena that affect the impact debris; however, such 
models are useful for several reasons. First of all, the 
models can be used to predict the temporal variation 
of the initially observed species; the degree to which 
the predictions match the observations tells us about 
the relative effectiveness of dynamical vs. photochem- 
ical processes in the post-SL9 Jovian stratosphere and 
can greatly enhance our understanding of complex pho- 
tochemical mechanisms that are unstudied in labora- 
tories on Earth. Secondly, the models can be used to 
direct the search for undetected molecular species and 
to identify the short- and long-term reservoirs of various 
elements. Finally, such models can aid in the analysis 
of the haze observations by predicting the composition, 
location, formation rate, and lifetime of possible con- 
densates in Jupiter’s stratosphere. 

Photochemical Model 

The hydrocarbon reaction list, absorption cross sec- 
tions, quantum yields, diffusion coefficients, boundary 
conditions, and background atmospheric structure of 
Gladstone et al. [1995] form the basis for our photo- 
chemical model, but we update the hydrocarbon pho- 
tochemistry and include oxygen, nitrogen, and sulfur re- 
actions. Our photochemical model contains 144 species 
and over 900 reactions. We have tried to identify and 
include the important reaction schemes that are most 
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likely to affect the abundances of the "‘parent” molecules 
observed after the SL9 impacts. However, the rates for 
several important reactions are not, available in the lit- 
erature, and we often had to make estimates. There- 
fore, this work represents a preliminary, qualitative at- 
tempt to follow the evolution of the important chemical 
species; the results should not be regarded as definitive. 
Our reaction list is available upon request. 

Since some of the observed species form most easily 
in low pressure re-entry shocks [Zahnle et al ., 1995], we 
assumed that most of the new material was deposited 
above 1 mbar in the stratosphere. Based on the results 
presented by Noll et al. [1994], we start with initial 
column abundances of reduced sulfur species as shown 
in Table 1. The initial abundances of the nitrogen- and 
oxygen-containing compounds are described in Moses et 
al. [1995]. Due to a lack of coherent information at this 
time, we do not enhance the hydrocarbon abundances 
after the plume splashdown, but begin with the back- 
ground hydrocarbon abundances derived by Gladstone 
et al. [1995] for the same latit ude region. The new post- 
impact photochemistry does not measurably affect the 
profiles of the major hydrocarbon compounds. 

Results 


The abundances of our initial sulfur species change 
with time as photochemical reactions are initiated. Ta- 
ble 1 illustrates the temporal variation of several of the 
more interesting species. Due to the large uncertainties 
in initial compositions at the impact sites, one should 
not put great reliance on the actual numbers in this ta- 
ble; two significant digits are included solely as an aid 
in comparing relative changes. 

Sulfur Photochemistry 


The most abundant sulfur compound observed im- 
mediately following the impacts is diatomic sulfur (S2)* 
S2 is quickly lost by photolysis in the upper atmosphere 
and by the formation of larger sulfur molecules in the 
lower stratosphere. In fact, the production of molecu- 
lar sulfur (Ss) by various three-body reactions is very 
rapid. We adopt the following scheme for the produc- 
tion of various sulfur molecules: 


s x A s x _, + s 
2S A S 2 
s + s 2 A s 3 

25 3 A S 6 
s 2 + s x A s x+2 

s + S x S 2 + S x _! 

254 A s 8 


for x = 2, 3, 4 


for x = 2, 4, 6 
for x = 3, 4, 5, 6 


where M represents any third molecule or atom. As Ta- 
ble 1 indicates, over 70% of the initial column budget of 
S2 molecules has been converted to Ss in the first hour 
after we begin our calculations, and another ~-T0% is 
converted to other sulfur molecules (e.g., S3, S4). Ss, 
which has a puckered ring structure, will condense al- 
most as soon as it is formed. 


On the other hand, Ss formation is only rampant in 
the middle and lower stratosphere at pressures greater 
than ~0.1 mbar. At higher altitudes, S2 is effectively 
photolyzed before larger sulfur molecules can form. The 
S2 (B-X) band system lies at 240 - 360 nm and is clearly 
observed by the HST FOS [Noll et al 1995]. At wave- 
lengths shorter than ~278 nm, S2 dissociates into two 
ground-state ( 3 P) sulfur atoms. The lifetime against 
photolysis for the S2 molecule at 10“ 5 mbar and 44° S 
latitude on Jupiter in our diurnally averaged model is 
5.3 hours. At high altitudes, photolysis seems to provide 
an effective permanent loss mechanism for the S2; how- 
ever, in the middle stratosphere near 0.01 - 0.1 mbar, 
the three-body reaction 2S+ M — »• S2 + M balances the 
photolysis, and S2 is efficiently recycled. Thus, we find 
that some S2 remains several months after the impacts 
(see Fig. 1 and Table 1). Because this result depends on 
the rates of several reactions that have not been mea- 
sured in the laboratory, this conclusion should be re- 
garded as speculative. The HST spectra show no evi- 
dence for S2 at the G site 22 days after the impact, but 
no upper limits are given in Noll et al. [1995], and it 
is unclear whether models such as are shown in Fig. 1, 
in which the S2 is confined to a narrow 7 altitude region, 
are consistent with these observations. 

Carbon disulfide (CS2), which is also observed un- 
ambiguously in the 1 1ST FOS spectra, has strong ab- 
sorption bands in the 185 - 230 nm region. At these 
wavelengths, CS2 dissociates into CS + S, where ~80% 
of the sulfur atoms emerge in the l D excited state and 
the rest in the 3 P ground state. The S( 1 D) reacts ef- 
ficiently with H2 to form SH + H. CS2, like S2, is not 
lost from the middle stratosphere as fast as its photol- 
ysis lifetime (8.2 hours at 10 -5 mbar) would indicate. 
Instead, CSo is partly recycled and is also formed dur- 
ing the photolysis of S2 in reaction schemes that involve 
S atoms and hydrocarbon radicals and have HvCS and 
HCS as intermediates. 

Thioformaldehyde (H2CS) forms in potentially ob- 
servable quantities if we assume (as we did) that the 



SULFUR (S 2 ) MIXING RATIO 


Figure 1. The evolution of S2. The solid line shows 
our assumed initial abundance, the dashed line shows 
the mixing ratio after 9 hours, the dot-dashed line after 
3 days, and the dotted line after 30 days. 
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unmeasured reaction S + CH3 — > H 2 CS + H is fast. It 
is possible, however, that S 4- CH3 will act to prefer- 
entially form CH3S via three-body reactions and that 
H2CS will not be a major product. If this alternative is 
true, then H 2 CS will cease to be a major component in 
our model while methyl mercaptan (CH3SH), dimethyl 
sulfide (CH3SCH3), dimethyl disulfide (CH3S) 2 , and or- 
ganic sulfur ring species (e.g., C2H4S, C3H6S) will be- 
come more important. CS production would also de- 
crease in this case. 

Carbon sulfide (CS) is produced by CS 2 photolysis 
and, more importantly in our model, by several different 
reaction schemes whose ultimate source is S 2 photolysis. 
Although CS is considered a radical species, very few re- 
actions seem to be effective at permanently removing it 
from Jupiter’s stratosphere. One exception is a reaction 
that we have invented: NH 2 + CS — > HNCS + H. We 
have assumed that this reaction proceeds relatively inef- 
ficiently, but the CS and the NH 2 radicals (which result 
from NH 3 photolysis) are so abundant that the reac- 
tion becomes important. As can be seen from Table 1, 
the abundance of isothiocyanic acid (HNCS) increases 
dramatically with time. We have not included any loss 
mechanisms for HNCS in our model, so if it does not 
condense in the stratosphere, it may be photolyzed be- 
fore significant amounts can build up. 

Our prediction for the CS column abundance 3 days 
after the impact is 2 x 10 15 cm“ 2 , which is an order of 
magnitude higher than the abundance derived by Noll 
et ai [1995] for CS observed in emission by the HST 
FOS three days after the G impact. However, our re- 
sults do compare well with the millimeter observations 
of Lellouch et ai [1995] also taken 3 days after the G 
impact if we assume that the ~3 x 10 11 g CS seen by 
these authors is spread over an areal spot ~7500 km in 
radius. The column quoted by Noll et ai would require 
the CS to be spread over a spot 25,000 km in radius in 
order to account for the Lellouch et ai observations. 

Hydrogen sulfide (H 2 S) is inferred, but not unam- 
biguously detected, in the HST FOS absorption spectra 
[ Yelle and McGrath , 1995; Noll et ai , 1995]. H 2 S has a 
nearly featureless absorption in the 165 - 250 nm region, 
with a broad peak located near 187 nm. Photolysis of 
H 2 S proceeds predominantly by H 2 S + hv -4 H -f SH; 
the HoS lifetime against photolysis is 1.8 days at 10 -5 


mbar. The SH radicals formed during the H 2 S pho- 
tolysis react with atomic H in the upper atmosphere 
to form S + H 2 and with various sulfur radicals in the 
lower stratosphere to form HS 2 . H 2 S is not efficiently 
recycled and is lost rapidly from the stratosphere. 

A close examination of Table 1 will reveal that ni- 
trogen-sulfur species are important reservoirs for the 
sulfur. This result is one of the most speculative of our 
paper. Sulfur-nitrogen reactions have not been exten- 
sively studied in laboratories, but both of the main N-S 
species in our model, NS radicals and HNCS molecules, 
have been observed in interstellar space. In our model, 
HNS radicals are produced through reactions of atomic 
sulfur with NH 2 . The HNS radicals react with both S 
and SH to form NS, but the only permanent loss for NS 
in our model is the fast reaction N + NS -A N 2 + S. Since 
N atoms are not abundant, this reaction is ineffective. 
As mentioned previously, HNCS is formed by the reac- 
tion of NII 2 with CS. If the reactions producing these 
species are ineffective or if these species arc readily lost 
by mechanisms we have omitted, then our conclusions 
regarding the major reservoirs of both sulfur and nitro- 
gen may change. The sulfur will still mostly end up in 
Sg, but this process may be accelerated, and ammonia 
(NH 3 ) may be more efficiently recycled and not lost so 
rapidly from our model [see Moses et a/., 1995]. 

Condensates 

Dark aerosol particles were immediately visible at the 
impact sites [e.g., West et ai , 1994]. These particles 
probably consist of silicates, metal sulfides or oxides, 
and/or organic debris (e.g., soot) that formed thermo- 
chemically in the impact or re-entry shocks. Does pho- 
tochemistry contribute to the aerosol properties? The 
particles themselves are very dark [ West et ai , 1995], 
and none of the species in our model form condensates 
that are as dark. In an interesting and relevant lab- 
oratory experiment, Khare et ai [1978] irradiated a 
mixture of CH 4 , NH 3 , and H 2 S gases with high doses 
of 254 nm ultraviolet radiation. The brown solid that 
formed during this procedure was composed of 84% Ss 
and 16% more complex organic matter. This experi- 
ment, and others like it, reveal that solids that contain 
a large percentage of Ss do not necessarily have to be 
bright at visible wavelengths. 


Table 1 . Time Variation of Column Abundances (cm 2 ) 


Species Initial 


1 hr 3 days 9 days 30 days 17 months 


S 2 

1.5 

x 10 18 

s* 


_ 

H 2 S 

7.7 

x 10 16 

h 2 s 4 


- 

CS 

1.4 

x 10 12 

cs 2 

2.0 

x 10 15 

h 2 cs 


- 

CH 3 SH 


- 

C 2 H 4 S 


- 

NS 


- 

HNCS 


- 


2.4 x 10 17 

2.7 x 10 17 

7.4 x 10 16 

6.0 x 10 13 

1.6 x 10 14 
1.9 x 10 15 

1.0 X 10 14 

2.1 x 10 u 

2.2 x 10 12 

2.3 x 10 14 

4.5 x 10 13 


1.2 x 10 17 

3.4 x 10 17 

4.4 x 10 15 
5.9 x 10 14 

2.3 x 10 15 

4.7 x 10 14 

1.2 x 10 15 

6.6 x 10 12 

3.8 x 10 12 

3.0 x 10 16 

7.7 x 10 15 


8.8 x 10 lb 

3.5 x 10 17 

2.0 x 10 14 

5.9 x 10 14 

5.7 x 10 15 

1.3 x 10 14 

9.9 x 10 14 

9.9 x 10 12 

3.3 x 10 12 

5.9 x 10 16 

2.0 x 10 16 


6.7 x 10 16 

3.5 x 10 17 

3.9 x 10 13 

5.9 x 10 14 

1.7 x 10 16 

4.1 x 10 13 

6.8 x 10 14 

1.5 x 10 13 

2.1 x 10 12 

6.8 x 10 16 

2.7 x 10 16 


8.9 x 10 15 

3.5 x 10 17 

6.2 x 10 13 

5.9 x 10 14 

1.1 x 10 17 

1.6 x 10 14 

1.0 x 10 14 

9.0 x 10 11 

1.2 x 10 11 

6.1 x 10 16 

3.2 x 10 16 


Times refer to approximate times after splashdown. Ss and H 2 S 4 will be primarily in the condensed phase. 
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West et al. [1995] claim that the visual properties 
of the aerosols do not change appreciably with time, 
except for changes that are consistent with the parti- 
cle size increasing due to coagulation. This result is 
somewhat surprising since one might naively expect the 
continual accumulation of Sg to brighten the aerosols 
over time — our predicted condensation flux for Ss (as 
much as 1 x 10 " 12 g cm -2 s” 1 ) is enormous. Perhaps 
small amounts of impurities have darkened the elemen- 
tal sulfur. Alternatively, signatures for Sg condensing 
about a darker particle might not show up clearly - 
pure Sg is very transparent at yellow-to-red and near- 
infrared wavelengths. Detailed modeling of the optical 
properties of relevant composite mixtures and a further 
examination of all the available data might help resolve 
this issue. Note that Sg absorbs strongly in the ultravi- 
olet; future photochemicals models should consider the 
effects of UV-shielding by the aerosols. 

Because we have assumed that the impact debris 
was mixed down to a few mbar at most and because 
the low eddy diffusion coefficients we adopt do not al- 
low for appreciable mixing in the lower stratosphere, 
the only molecules that can condense in our model are 
species of low volatility and reasonable abundance. The 
only sulfur-containing molecules that, fit this descrip- 
tion (other than Ss) are H 2 S 4 and higher-order sul- 
fanes. We were unable to discover a molecule that forms 
dark aerosols yet that is volatile enough to account for 
the transient atmospheric ring features that several re- 
searchers suggest might be due to condensation and 
evaporation during the passage of an atmospheric w r ave. 

Summary and Conclusions 

In our model, sulfur photochemistry at the SL9 im- 
pact sites is rapid and complex. Condensed Sg is the 
major reservoir for sulfur after the first, few hours, but 
CS and nitrogen-sulfur species may also become impor- 
tant after several days, and the inventories of sulfanes 
(H 2 s x ) and thioformaldehyde (H 2 CS) are nontrivial. 
The observed condensates should be rich in sulfur. The 
photochemical model predictions from Table 1 can be 
used in conjunction with observations to shed light on 
important photochemical processes in the post-impact 
Jovian stratosphere. For instance, we predict that S 2 
will persist in the stratosphere longer than its photol- 
ysis lifetime warrants. HST FOS spectra [Noll el al , 
1995] indicate the presence of S 2 three days after im- 
pact but not 23 days later. Is our narrowly-confined S 2 
profile consistent with these observations? If not, can 
that result tell us something about the rates of some un- 
measured sulfur reactions? The CS variations observed 
at millimeter wavelengths [Lellouch et a/., 1995] are af- 
fected by horizontal spreading as well as absolute abun- 
dances. Our results might help in the analysis of this 
complex behavior. Our prediction that CS 2 is lost more 
rapidly than NH 3 in the first month after the impacts 
[see Moses et a/., 1995] is consistent with HST FOS ob- 
servations [Yelle and McGrath , 1995]; however, we may 
be losing CS 2 more rapidly in our model than it is lost 


on Jupiter. All of our conclusions are highly specula- 
tive and depend on several kinetic reaction rates and 
initial conditions that are unknown at this time. We 
will examine the sensitivity of our model to different 
parameters in a later work. 
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